Distorting the Hump-backed Plateau of Jets with Dense QCD Matter 
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The hump-backed plateau of the single inclusive distribution of hadrons inside a jet provides a 
standard test of the interplay between probabilistic parton splitting and quantum coherence in QCD. 
The medium-induced modification of this QCD radiation physics is expected to give access to the 
properties of the dense medium produced in relativistic heavy ion collisions. Here, we introduce a 
formulation of medium-induced parton energy loss, which treats all leading and subleading parton 
branchings equally, and which - for showering in the vacuum - accounts for the observed distribution 
of soft jet fragments. We show that the strong suppression of single inclusive hadron spectra 
measured in Au-Au collisions at the Relativistic Heavy Ion Collider (RHIC) implies a characteristic 
distortion of the hump-backed plateau; we determine, as a function of jet energy, to what extent the 
soft jet fragments can be measured above some momentum cut. Our study further indicates that 
the approximate p-r-independence of the measured nuclear modification factor does not exclude a 
significant Q 2 -dependence of parton energy loss. 



Over the last five years, experiments at RHIC have 
established a phenomenon of strong high-p^ hadron sup- 
pression [l|, |2|, which supports the picture that high- 
Pt partons produced in the dense matter of a nuclear 
collision suffer a significant energy degradation prior to 
hadronization in the vacuum Q. The microscopic dy- 
namics conjectured to underly high-pr hadron suppres- 
sion is medium-induced gluon radiation 0,0, ||J, i.e., a 
characteristic medium- induced distortion of the standard 
QCD radiation pattern tested extensively by jet mea- 
surements in high energy e + e~ and pp (pp) collisions. 
Modeling this effect accounts for the measured single- 
inclusive spectra and leading back-to-back hadron corre- 
lations However, existing models indicate only 
rough qualitative features of subleading jet fragments 
such as their broadening and softening due to medium- 
effects. This is so mainly, since medium-induced parton 
splitting is included for the leading partons only, but not 
consistently for the subleading splitting processes 0, , 
and since energy momentum conservation is taken into 
account globally, but not ensured locally for each parton 
splitting 0, 0. Here, we propose a formulation which 
overcomes these limitations. 

Subleading jet fragments are known to provide many 
fundamental tests of QCD radiation physics. In partic- 
ular, for soft particle momentum fractions x — p/Ej Ct 
inside a quark- or gluon- [i = q, g) initiated jet of en- 
ergy and virtuality Q ~ -Ejet, the single inclusive distri- 
bution Di(x,Q 2 ) is dominated by multiparton destruc- 
tive interference, and thus tests quantitatively the un- 
derstanding of QCD coherence |10l lllj . Remarkably, to 
double and single logarithmic accuracy in £ = In [1/ x] 
and t = ln[Q/A e ff], A c g — O(Aqcd), the effects of 
this destructive quantum interference can be accounted 
for by an angular ordering prescription of a proba- 
bilistic parton cascade with leading order (LO) split- 
ting functions. The so-called modified leading logarith- 



mic approximation (MLLA) resums these effects and ac- 
counts for the lar ge ^/ aj next-to-leading corrections of 
Di{x,Q 2 ) HI [Hill. The MLLA leads to an evolu- 
tion equation for the v-th Mellin moments Mi(v,r) = 

/ °°^ e -^xA(x,Q 2 ) [HEEl, 
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Here, $jj denote combinations of particular moments of 
leading order splitting functions, for example 
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$,,(!/) = 2 / dzP qq (z) (z v -1) 



The shift {y + ■§-) in (JIJ accounts for angular ordering. 
For a parton fragmentation which starts at high initial 
scale t and ends at some hadronic scale To, the solution of 
JTJ has to fulfill the initial conditions M(x, tq) = 6(1 — x) 
and JppM(x,r — r ) = 0, since the parton must not 
evolve if produced at the hadronic scale. 

The lowest Mellin moments v ~ determine the main 
characteristics of Di(x,Q 2 ). For an approximate solu- 
tion of l|T)l. one can thus expand the matrix in Q to 
next-to-leading order in (y + J^) and diagonalize it. Its 
eigenvalue with leading 1/ [y + J^)-term yields a differ- 
ential equation of the confluent hypergeometric type ^lj • 
This leads to an analytic expression for D(x,Q 2 ), whose 
shape does not distinguish between quark and gluon par- 
ents, since the multiplicity is dominated in both cases by 
gluon branching. For the hadronic multiplicity distribu- 
tion dN h /d^, one assumes that at the scale To, a parton is 
mapped locally onto a hadron with proportionality factor 
K h ~ 0(1) ("local parton hadron duality", LPHD) 
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Com par isons of to data have been performed repeat- 
edly |l3l llil llol Ha . Il^l over a logarithmically wide kine- 
matic regime 7 < Ej et < 150 GeV in both e + e" and 
pp/pp collisions. To illustrate the degree o f ag reement, 
we reproduce in Fig. Q two sets of data 0, Il6j together 
with the curves obtained from ©. The parameters K h 
and A c ff entering (J3J) were chosen as in Refs. 0, llij . 
A cfr = 254 MeV, K h = 1.15 for E jct = 100 GeV, 
K h = 1.46 for Ejet = 7 GeV. Following Ref. 0, we 
use iV/ = 3. From Fig. we conclude that Eq.© ac- 
counts reasonably well for the jet multiplicity distribution 
in the kinematic range accessible in heavy ion collisions 
at RHIC (Ejet ~ 10 GeV) and at the LHC (E jct ~ 100 
GeV). Corrections not included in © are of relative or- 
der 1/r, which at face value corresponds to a 30% (15%) 
uncertainty at typical RHIC (LHC) jet energies. Also, 
the MLLA resums large £, r ~ £, but is expected to be 
less accurate for hard jet fragments, where other improve- 
ments are currently sought for [lsj . Thus, the agreement 
of to data for the entire £-range is surprisingly good. 
At least from a pragmatic point of view, © can serve 
as a baseline on top of which one can search for medium 
effects. 
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the gluon bremsstrahlung are more singular than <i/ vac 
for small us and may thus be expected to dominate the 
multiplicity distribution J3J). However, destructive inter- 
ference due to finite in-medium path length is known to 
regulate the soft w-divergence (2TJ . For the relevant range 
of soft us, this may be modeled as us 3 ^ — ~ / mo d = const. 



A medium-induced gluon bremsstrahlung spectrum, con- 
sistent with this ansatz, was also found in 21]. This 
suggests that medium effects enter (01 by enhancing the 
singular parts of all LO splitting functions P gg , P qg , P qq 
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We do not modify the non-singular subleading terms. 
On general grounds, one expects that medium-induced 
rescattering is a nuclear enhanced higher-twist contribu- 
tion (/med ~ gk) H3- This means that it is subleading 
in an expansion in Q 2 , while being enhanced compared 
to other higher twist contributions by a factor propor- 
tional to the geometrical extension ~ L of the target. A 
l/Q 2 -dependence of / me d is also suggested by the follow- 
ing heuristic argument. A hard parton of virtuality Q has 
a lifetime ~ 1/Q in its own rest frame, and thus a life- 
time (in- medium path length) t = ^ q before it branches 
in the rest frame of the dense matter through which it 
propagates. Medium effects on a parton in between two 
branching processes should grow proportional to (some 
power of) the in-medium path length and thus oc l/Q 2 
or higher powers thereof. 

In contrast, jet quenching models 0,0,013 reproduce 
inclusive hadron spectra in Au-Au collisions at RHIC by 
supplementing the standard QCD LO factorized formal- 
ism with the probability P(AE) that the produced par- 
tons radiate an energy AE due to medium effects prior 
to hadronization in the vacuum |9|| 



FIG. 1: The single inclusive hadron distribution as a function 
of £ = In [j?jet /p] . Data taken from e + e~ collision experiments 
TASSO [13 and OPAL p]J, E iet = vV 2 - Lines through data 
obtained from the MLLA result ||3J . Dashed and dash-dotted 
curves labeled "in medium" are calculated with a medium- 
modification / mc d = 0.8 of the LO splitting functions. 



The multiplicity distribution dN h /dt; is dominated by 
soft gluon bremsstrahlung, dl vac ~ C R ,^^^k 
lu = z Ej ct , which is described by the singular parts 
~ |, ~ (i—z) °^ ^ e Q^D splitting functions entering 
(J2J . They determine the leading i-terms of the evolution 
matrix in QJ. Remarkably, calculations of the additional 
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medium-induced radiation indicate that 
if the medium is modeled by soft multiple momentum 
transfers 0,H3], and ~ ^ if the medium is modeled by a 
single hard momentum transfer 0, l2fj| . Thus, parametri- 
cally, the additional medium-dependent contributions to 



OO " U r. 



dl med (uii) 



n=0 Li=l 

- poo d j med 

x exp — / dus 

L Jo dus 



dus 



s[aE-^ 



(5) 



This formula is based on a probabilistic iteration of 
medium-modified parton splittings, but does not keep 
track of virtuality or angular ordering. The fey-integrated 
medium-induced contribution dl mcd is treated on an 
equal footing with LO vacuum splitting functions. In 
this sense, the medium-modified fragmentation func- 
tion D { ™ d \ Xl Q 2 ) = Jl deEP(AE) ^ D h/ Jp,Q 2 ), 
e = AE/E, entering jet quenching models [HEl 0> El> 
amounts to a medium-induced Q 2 -independent modifi- 
cation of parton fragmentation. 

The single inclusive distribution D(x,Q 2 ), supple- 
mented by LPHD, is a fragmentation function. Single 
inclusive hadron spectra, whose parent partons show a 
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power law spectrum oc l/p T ^ PT \ test D(x,Q 2 ) in the 
range, in which x nt ^ PT ^ 2 D(x,Q 2 ) has significant sup- 
port. However, for large x, the accuracy of the MLLA 
result for D(x,Q 2 ) becomes questionable 0, To 
understand to what extent the MLLA result may still be 
used from a practical point of view, we have compared it 
to the KKP parametrization of fragmentation func- 
tions. In the range of Q 2 and x relevant for single inclu- 
sive spectra (0.4 < x < 0.9), we observe that the KKP 
and MLLA fragmentation functions both drop by 2 or- 
ders of magnitude. They do show somewhat different 
shapes but - after adjusting the overall normalization - 
they differ for all a;-values by maximally ~ 30% or signif- 
icantly less (data not shown). For the nuclear modifica- 
tion factor Raa, which is the ratio of modified and un- 
modified single inclusive hadron spectra, and which does 
not depend on the overall normalization oi D(x, Q 2 ), this 
is a relatively small uncertainty, if one aims at character- 
izing a factor 5 suppression. We thus conclude that the 
MLLA fragmentation function obtained from (JJJ can be 
used to calculate Raa- 

Saa(Pt) Au-Aih;t , 0-10% centrality 
1 
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FIG. 2: The pT-dependence of the nuclear modification factor 
Raa, calculated for a medium-enhanced parton splitting with 
/mod- Data taken from [2^1 . 



To determine Raa, we have parametrized the par- 
tonic pr-spectrum at RHIC energies ^/smn = 200 GeV 
by a power law l/p T (pT \ n{p T ) = 7 + 0.003p|/GeV 2 , 
which accounts for kinematic boundary effects at pj< ~ 
0(-Jsnn). Single inclusive hadron spectra and Raa are 
calculated by convoluting this spectrum with D(x,Q 2 ). 
Medium effects are included through the factor / me d in 
the singular parts of all LO parton splitting functions, see 
Eq.lUJ. As seen in Fig.0 the choice / mc d — 0.6 -=- 0.8 re- 
produces the size of the suppression of Raa ~ 0.2 in cen- 
tral Au-Au collisions at RHIC |2^| . Jet quenching models 
based on JSJ) yield a slightly increasing pr-dependence of 
Raa{pt) for a power law n(pr) = n, and a rather flat de- 
pendence if trigger bias effects due to the pr-dependence 
of u{pt) are included I n contrast, the MLLA re- 

sult for Raa{pt) decreases with increasing px, and this 



tendency is even more pronounced for a realistic shape 
of the underlying partonic spectrum, see Fig. [2] The 
reason is that in the MLLA, parent partons of higher 
Pt have higher initial virtuality Q ~ px, and undergo 
more medium- induced splittings; this results in a smaller 
value of Raa- A proper treatment of nuclear geometry 
may affect quantitative aspects of Fig. but is unlikely 
to change this qualitative observation. Hence, the ob- 
served flat pr-dependence of Raa (jpt ) , one may require 
(in accordance with the arguments given above) a non- 
vanishing Q 2 -dependence of / mc d, which would reduce 
medium-effects on high-py (pr ~ Q) partons. 

Motivated by this observation, we have attempted 
to solve Eq.lQJ for a non-trivial Q 2 -dependence of the 
medium-enhancement / mc d- We did not find an analyti- 
cal solution. However, in the absence of medium effects, 
the analytical solution J3J is reproduced by Monte Carlo 
(MC) parton showers based on angular ordering fl3| . 
This remains true for non- vanishing / mo d- The present 
study can serve to check future MC showers implement- 
ing Q), and it can be extended in MC studies to include 
a non-trivial Q 2 -dependence of /med- We plan such a 
MC study, mainly to establish to what extent the ap- 
proximate pT-independence of Raa up to pt ~ 10 GeV 
allows for a significant Q 2 -dependence of parton energy 
loss. The question of whether and on what scale these 
effects are l/Q 2 -suppressed is of obvious importance for 
heavy ion collisions at the LHC, where medium- modified 
parton fragmentation can be tested in a logarithmically 
wide Q 2 -range. 

What is the distortion of the longitudinal jet multi- 
plicity distribution @, consistent with the observed fac- 
tor ~ 5 suppression of Raa? In contrast to calcula- 
tions based on (|5J), the medium-enhanced parton split- 
ting introduced via MLLA conserves energy-momentum 
exactly at each branching, it treats all secondary branch- 
ings of softer gluons equally, and it continues all branch- 
ings down to the same hadronic scale. This makes it 
a qualitatively improved tool for the calculation of lon- 
gitudinal multiplicity distributions, since it matters ob- 
viously for dN/dt; whether one gluon is radiated into 
the bin £ = In [-Ejct/Pgiuon], or - after further splitting 
g — > g(z) g(l — z) — two gluons with momentum fractions 
z and (1—z) into bins £+ln [1/(1 — z)], £+ln [1/z], respec- 
tively. We have calculated dN/d£ for a medium-enhanced 
parton splitting / me d = 0.8 consistent with Raa = 0.2. 
Results for jet energies relevant at RHIC and at the LHC 
are shown in Fig. ^ In general, the multiplicity at large 
momentum fractions (small £) is reduced and the cor- 
responding energy is redistributed into the soft part of 
the distribution. The maximum of the multiplicity dis- 
tribution also shifts to a softer value, but this shift is 
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Many experimental characterizations of the medium- 
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within jets, which have been studied in the absence of 
medium effects IT2I ll3T ] . It may also apply to transverse 
jet broadening |l2j, which for a Q 2 -dependent / mo d may 
be significantly reduced since large angle emission would 
remain essentially unmodified by the medium. We plan 
to address these open questions in future work. 

We thank S. Catani, G. Sterman and B. Webber for 
helpful discussions. 



FIG. 3: The change of the total hadronic multiplicity inside 
jets of different energy -Ejct, as a function of the soft back- 
ground cut p^ 1 * above which this multiplicity is measured. 
Medium-effects are modeled by / mc< j = 0.8. 

modified internal jet structure in heavy ion collisions at 
RHIC and at the LHC require a soft momentum cut pip* 
to control effects of the high multiplicity background. 
Can one observe the increase in soft multiplicity shown in 
Fig. if such a soft background cut p!p* is applied? To 
address this issue, we have calculated from © the total 
hadronic multiplicity N h (pt > p™') above p™*. As seen 
in Fig. the medium-enhanced component of soft mul- 
tiplicity lies below a critical transverse momentum cut 
Profit which increases significantly with E- ]ct . For a typ- 
ical hard jet at RHIC (£j e t = 15 GeV), the additional 
soft jet multiplicity lies buried in the soft background, 
p^ rit ~ 1.5 GeV. For E jct = 100 GeV, accessible at the 
LHC, Pr U crit — 4 GeV lies well above a cut which de- 
pletes the background multiplicity by a factor 10. For 
E iet = 200 GeV, we find p§?* rit ~ 7 GeV. The associated 
total jet multiplicity N h (pr > P^ci-it) f° r these jet ener- 
gies rises with £j et from ~ 4, to ~ 7. Fig. [3] indicates a 
qualitative advantage in extending jet measurements in 
an LHC heavy ion run near design luminosity to signifi- 
cantly above E- iCt ~ 100 GeV, where a sizeable kinematic 
range 2-^3 GeV < pt < Px^rit — ^ GeV becomes ac- 
cessible. This may allow a detailed characterization of 
the enhanced medium-induced radiation above the soft 
background. 

In general, the formulation of parton energy loss within 
the MLLA formalism allows one to address several funda- 
mental questions, that remain untouched by recent model 
studies of jet quenching. This concerns in particular the 
important issue of the Q 2 -dependence of parton energy 
loss discussed above. Moreover, the use of a probabilistic 
formulation based on angular ordering can also be viewed 
as a test of the unproven assumption that the medium- 
induced destructive interference of multi-parton emission 
can indeed be accounted for by angular ordering, in close 
similarity to gluon radiation in the vacuum. We finally 
note that the formalism introduced here is not limited 
to a discussion of the hump-backed plateau: e.g. it can 
be applied to the calculation of two-particle correlations 
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